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Ion channels Background: Biological neuron models mainly analyze the behavior of neural networks. Neurons are described 

Postsynaptic neuron in terms of firing rates viz an analog signal. Purpose: The Izhikevich neuron model is an efficient, powerful 

synapse model of spiking neuron. This model is a reduction of Hodgkin-Huxley model to a two variable system and 

synaptic cleft is capable of producing rich firing patterns for many biological neurons. Methods: In this paper, the Regular 

postsynaptic membrane Spiking (RS) neuron firing pattern is used to simulate the spiking of Glutamate gated postsynaptic membrane. 

Glutamate Results: Simulation is done in MATLAB environment for excitatory action of synapses. Conclusions: Analogous 

ENFET simulation of spiking of excitatory postsynaptic membrane potential is obtained. 
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Introduction 

A typical neuron can be divided into three functionally distinct 
parts: dendrites, soma, and axon. Dendrites play the role of 
the 'input device' that collects signals from other neurons and 
transmits them to the soma. The soma is the central processing 
unit. If the total input exceeds a certain threshold, then an out- 
put signal is generated. The junction between two neurons is 
called a synapse. The neuronal signals consist of short electrical 
pulses. The pulses, so-called action potentials or spikes, have an 
amplitude of about 100 mv and typically a duration of 1-2 ms. 
The form of the pulse does not change as the action potential 
propagates along the axon. A chain of action potentials emit- 
ted by a single neuron is called a spike train - a sequence of 
stereotyped events which occur at regular or irregular intervals. 
Since all spikes of a given neuron look alike, the form of the ac- 
tion potential does not carry any information. Rather, it is the 
number and the timing of spikes which matter. Action poten- 
tials in a spike train are usually well separated. Even with very 
strong input, it is impossible to excite a second spike during or 
immediately after a first one. 1 

One of the most important neuron models in computational 
neuroscience is the Hodgkin-Huxley model. Based on their 
1 952 experiments on the giant axon of the squid, Hodgkin and 
Huxley developed a mathematical model describing neuron 
membrane potential and the flow of ions through channels 
in the membrane. In the Hodgkin-Huxley model (Figure 1) the 
cell membrane is considered as a capacitor. An input current 
I applied to the membrane may be split into a current which 
charges the capacitor l cap and a current that passes through the 
ion channels I u } 



I : 



cap chan 



The ionic current is divided into components carried by sodium 
and potassium ions l Na and l K respectively, and a small 'leakage 
current' (l 0 ) made up by chloride and other ions. Each compo- 
nent of the ionic current is determined by a driving force which 
may conveniently be measured as an electrical potential differ- 
ence and a permeability coefficient which has the dimensions 
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Fig. 1: H-H model. 

of a conductance. Thus the sodium current (l Na ) is equal to the 
sodium conductance (g Na ) multiplied by the difference between 
the membranes potential (E) and the equilibrium potential for 
the sodium ion (E Na ). The experiments suggest that g Na and gk 
are functions of time and membrane potential, but that E Na , 
E k , E 0 , C M and g o may be taken as constant. 3,4 
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=C M (dV m /dt) + g o (V m -E o )-g Na (V m -E Na ) + g k (V m -E k ) 

where V m represents the postsynaptic membrane potential es- 
tablished by the ionic and capacitive membrane current, C M is 
the capacitance of the lipid bilayer of postsynaptic membrane, 
t is time. The Hodgkin-Huxley model is biophysically meaning- 
ful but very expensive to implement. The motive of this paper is 
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to develop a simple analog circuit model that can simulate the 
spiking of glutamate gated postsynaptic membrane. 

Glutamate sensitive ENFET 

Glutamate is a non essential amino acid. It is the primary excit- 
atory neurotransmitter in the human central nervous system. 
In simplest case, the binding reaction may be represented as: 

K 

Glutamate + Receptor(closed) <-> Glutamate-Receptor(open) 

k 2 equation (1) 

where k 1 and k 2 are the forward and backward rate constants 
respectively. 5 The field effect transistor (FET) gate surface plays 
an important role in the sensitivity and stability of the sensor. 
Each surface layer possesses certain pH sensitivity and can, 
therefore, detect minute changes in pH close to the electro- 
lyte/insulator interface. Tantalum pentoxide (Ta 2 0 5 ) is a promis- 
ing gate oxide material for sensoric purposes, as it has a large 
number of surface sites that leads to a large buffer capacity. 6 
The glutamate sensitive ENFET is prepared by immobilizing glu- 
tamate oxidase on the surface of gate oxide (Ta 2 0 5 ) (Figure 2). 
It is based on the biocatalyzed hydrolysis of L-glutamate in the 
presence of glutamate oxidase in accordance with the chemical 
reaction: 




L«glLtamate 2-woglLrtanate 

The enzymatic reactions on the modified electrode surface in- 
volved in the detection of glutamate are as follows: 7 

L-Glutamate + 0 2 ^2-Oxoglutarate + H 2 0 2 

H 2 0 2 ^0 2 + 2H++2e- 

The proton generated in this reaction changes the pH inside the 
enzyme which is registered by the underlying ion sensitive FET. 
The threshold voltage of such device, V TH (IS), is a function of 
pH of solution dependent on the concentration of glutamate. 
For very small value of drain to source voltage of ENFET, V ds , the 
conductance of such ENFET can be expressed as: 

G ds = W V gs- v TH ( |S )) equation (2) 

(3 is the geometric sensitivity parameter given by 
f3 = ^xC 0X (W/L 

where C ox is the oxide capacity per unit area, W and L are the 
width and the length of the channel respectively, and a, is the elec- 
tron mobility in the channel. V gs is the voltage applied to the refer- 
ence electrode and V TH (IS) is the threshold voltage of the ENFET. 
V TH (IS) in equation (2) can, therefore, be modeled as: 89 

V TH (IS)(t) = V^JO-exp^t) + exp(-k 2 t)U(t-t m )] equationp) 
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Fig. 2: Glutamate ENFET (a) Schematic diagram (b) Electronic diagram. 

where k.,and k 2 are time constants analogous to the rate con- 
stants of equation (1 ), U(t-t m ) is the Heaviside function and V TH0 
is the threshold voltage proportional to the maximum attain- 
able conductance, when all the transmitter-gated channels for 
Na + ions are open. 

Spiking neuron model 

The simplest form of spiking neural model includes time and 
gives their output in the form of spikes. A spiking model is a 
mathematical model which describes how input spike trains 
(sequences of timings) are mapped to an output spike train. 
Thus the output can be characterized by 

S = (t:i = 1.2 n), ti< t + 1 

Where t is the i th spike train in a train of n spikes. 10 

In this paper, a simple spiking model represented by equation 
(4), (5) is as biologically plausible as the Hodgkin-Huxley model, 
yet as computationally efficient as the integrate-and-fire model. 
The model is a reduction of the Hodgkin-Huxley model to a two 
variable system. Two ordinary differential equations describe 
the membrane potential v and a recovery variable u by: 
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dv/dt = 0.04v 2 + 5v + 1 40-u + I equation (4) 

du/dt = a(bv-u) equation (5) 

with an after-spike reset rule: 

if v > 30 then c, u^ u + d equation (6) 

Here v and u are dimensionless variables, and a, b, c, and d are 
dimensionless parameters. The variable v represents the mem- 
brane potential of the neuron and u represents a membrane 
recovery variable, which accounts for the activation of K+ ion- 
ic currents and inactivation of Na+ ionic currents, and it pro- 
vides negative feedback to v. After the spike reaches its apex 
(+30 mV), the membrane voltage and the recovery variable 
are reset according to the equation (6). Synaptic currents or in- 
jected dc-currents are delivered via the variable I. 11 The model 
can exhibit firing patterns of all known types of cortical neu- 
rons with the choice of parameters: a, b, c and d. In this paper, 
regular spiking neuron is used to generate spike. The Regular 
spiking (RS) neuron (Figure 3) has four model parameters a, b, 
c and d. The parameter a describes the time scale of the recov- 
ery variable u. Smaller values result in slower recovery. Typical 
value of a is 0.02. The parameter b describes the sensitivity of 



regular spiking (RS) 




Fig. 3: Showing the firing pattern of a regular spiking neuron. 
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Fig. 4: Schematic diagram of a single spike together with Afterhyperpo- 
larization (AHP) and afterdepolarization (ADP). 



the recovery variable u to the subthreshold fluctuations of the 
membrane potential v. Its value is 0.2. The parameter c de- 
scribes the after-spike reset value of the membrane potential v 
caused by the fast high-threshold K+ conductances. Its value is 
c = -65 mV. The parameter d describes after-spike reset of the 
recovery variable u caused by slow high-threshold Na+ and 
K+ conductances. A typical value is d = 8. 12 

Methods 

RS neurons generate 'regular' action potentials. Their synaptic 
function is excitatory. The slow potential of the RS neuron i 
begin to rise due to the input from some other neurons. When 
the potential reaches at the threshold of the RS neuron, the 
neuron fires. Immediately after the firing, the slow potential 
\J. is lowered to a reset value U R which is lower than the rest 
potential U o and after passing the period of absolute refractori- 
ness T RS it resumes rising again. This process of lowering and 
rising of the potential represents afterhyperpolarization (AHP) 
and afterdepolarization (ADP) observed in the membrane po- 
tential of RS neuron (Figure 4). 13 

The strategy that is used to produce spiking in the glutamate 
gated postsynaptic membrane is governed by a) integrating set 
of differential equations with incoming spikes inducing dis- 
crete changes in the state variables and outgoing spikes are 
triggered by a threshold condition b) using Euler's forward ap- 
proximation method. The modeling for excitatory synapse is 
shown in Figure 5. The leakage current l Q is considered to be 
small enough to be neglected. Since only sodium channels are 
responsible for excitatory action, the postsynaptic membrane 
is divided into three patches to represent spatial summation of 
the sodium current controlled by 



I = C(dV„/dt) + g o (V m -E o )-g Na (V m -E Na ) + g K (V m -E K ) 

where g Na is the total sodium conductance and g K is the non- 
gated potassium conductance. Vg 1# Vg 2 and Vg 3 are the volt- 
ages applied to the reference electrodes of the ENFETs. The 
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Fig. 5: Circuit model for Postsynaptic membrane. 



membrane potential V m is obtained by spatially and temporally 
varying g Na of glutamate-gated sodium channels. 

Simulation 

The component values assigned in the model for MATLAB 
simulation are taken from reference: 14 C m = 1 uf per cm 2 , 
g K = 36 mS per cm 2 , E Na = 1 1 5 mV and E K = -12 mV and I 
0. The specifications for three n-channel ENFETs are L = 1 5 |im f 
W = 2 [i,m, t ox — 100 nm, \l — 600 cm 2 /V-sec. The parameters 
for exponential function in equation (3), applied to each ENFET 
inputs are: V TH0 = -10 mV, t m — 600 i^sec, k 1 = k 2 = 5 msec. The 
three gates to source voltage of three ENFETS i.e Vg 1# Vg 2 and 
Vg 3 are kept constants at 1 Volt each. The three input param- 
eters of ENFET namely V TH1 , V TH2 and V TH3 dependence on con- 
centration of glutamate are applied in a staggered sequence 
at 0.01 msec intervals. This is done to simulate the time varia- 
tion in glutamate transmitter -receptor binding with respect 
to different patches of postsynaptic membrane. The values as- 
signed for regular spiking neuron model parameters: a = 0.02, 
b = 0.2, c = -65, d = 8. 

Results 

The MATLAB simulation outputs are shown in figure 6. 
The waveform represents the normal postsynaptic membrane 
potential with respect to time. Membrane potential, V m is 
established by spatial summation and temporal integration 
of the glutamate-gated sodium current. In this model, ac- 
tion potential takes the form of spikes and occurs during the 
time period of the pulse. Here the action potential is ex- 
hibited whenever V m reaches the threshold in the range of 
-60 mv and after that the action potential is reset when it 
reaches 40 mv. 

Conclusions 

The figure above shows an analogous simulation showing the 
response of the model to a time varying injected current. The 
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Fig. 6: Simulated result of excitatory postsynaptic membrane potential. 



three glutamate gated sodium channels are staggered at 0.01 ms 
time interval which simulates the time variation in transmitter- 
receptor binding with respect to different patches of post synap- 
tic membrane. The work shows that glutamate-sensitive ENFET 
can be used as a circuit analog to simulate the spiking of excit- 
atory postsynaptic membrane potential. 
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